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n-log-n  algorithm 


facility  location 
efficient  locations 


Given  n planar  existing  facility  locations,  a planar  new  facility  location 
X Is  called  efficient  If  there  Is  no  other  location  T at  least  as  close  to 
every  existing  facility  as  X,  end  ei^rlctly^^oser  than  X to  at  least  one  exist- 
ing facility.  iinneul  an  algorl'th&^miSi  Is  cither  of  order  n(log  n)  or 
order  n (depending  upon  hoir  the  problea  Is  defined)  that  constructs  all  effi- 
cient locations,  end  establish  that  no  eltemetlve  algoritha  can  be^pf  a lower 
order.  With  the  exception  of  two  coaputatlonel  coaplexlty  results , work 


Is  entirely  self-contained,  and  relies  almost  entirely  upon  simple  geometrical 
analyses . 
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Given  n planar  existing  facility  locations,  a planar  new  facility  i 

location  X is  called  efficient  if  there  is  no  other  location  Y at  least  * 

as  close  to  every  existing  facility  as  X,  and  strictly  closer  than  X to  ^ 

at  least  one  existing  facility.  We  present  an  algorithm  which  is  either  | 

j 

of  order  n(log  n)  or  order  n (depending  upon  how  the  problem  is  defined) 
that  constructs  all  efficient  locations,  and  establish  that  no  alternative 
algorithm  can  be  of  a lower  order.  With  the  exception  of  two  computational 
complexity  results,  our  work  is  entirely  self-contained,  and  relies  almost 
entirely  upon  simple  geometrical  analyses. 


INTRODUCTION 


Suppose  a number  of  existing  facilities  are  given,  having  planar 
locations  P^,  P^.  A new  facility  is  to  be  located  in  the  plane 

at  some  point  X to  be  determined.  With  P^  - (a^^,  b^) , X ■ (x,  y)  we 
denote  the  rectilinear  distance  between  X and  P^  by  r(X,  P^)  where, 
by  definition. 

r(X,  Pj)  - |x  - a^ I + |y  - b^|. 

Given  any  point  Y in  the  plane  for  which 
r(Y,  P^)  £ r(X,  P^),  1 1 i 1 n, 

we  say  that  Y dominates  X if  X and  Y satisfy  the  n inequalities,  with 
at  least  one  inequality  holding  strictly.  In  other  words,  Y dominates 
X if  Y is  at  least  as  close  to  every  existing  facility  location  as  X is, 
and  closer  than  X to  at  least  one  existing  facility  location.  If  no 
point  in  the  plane  dominates  a point  X*,  we  say  that  X*  is  efficient. 

We  denote  the  set  of  all  efficient  points  by  S*,  and  call  S*  the  efficient 
set.  Note  that  each  existing  facility  location  Pj  is  in  S*,  for  to  have 
a point  Y dominate  P^  would  imply  r(Y,  P^)  r(Pj , P^)  - 0,  in  turn 
Implying  Y ■ P^ . 

Recently  Wendell,  Hurter,  and  Lowe  [8]  have  Introduced  and  studied 
the  problem  of  finding  S*,  and  discuss  some  application  contexts,  with 
emphasis  upon  multiple  objective  problems.  We  remark  that  S*  may  also  be 
of  value  in  carrying  out  sensitivity  analyses  for  single  objective  location 

problems,  since  such  problems  typically  have  the  property  that  their 

* 

optiMl  solutions  are  efficient.  Also  it  appears  that  S may  be  of  value 
in  the  study  of  some  internal  warehouse  location  problems,  in  which  case 
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each  existing  facility  would  be  a warehouse  dock,  and  the  new  facility 
location  would  be  the  location  of  an  item  in  the  warehouse. 

Wendell,  Hurter,  and  Lowe  develop  a nuadier  of  properties  of  S*,  and 
present  two  different  algorithms  for  constructing  S*.  Their  work  relies 
upon  a good  deal  of  relatively  deep  convexity  analysis  [2  ],  [7  ]•  We 
establish  in  this  paper  that  S'*  can  be  characterized  in  an  entirely  self- 
contained,  simple,  and  intuitively  appealing  manner,  using  only  geometry. 

We  consider  our  work  both  complements  and  supplements  the  work  of  Wendell, 
Hurter,  and  Lowe.  In  particular,  we  point  out  that  the  arrow  algorithm 
we  present  for  constructing  S*  is  closely  related  to,  and  motivated  by, 
the  second  of  the  two  algorithms  in  [ 8 ] • 

The  primary  value  of  the  arrow  algorithm  In  our  approach  is  as  a tool 
to  facilitate  proofs.  A second  algorithm  we  present,  the  row  algo- 
rithm, is  more  efficient.  In  fact,  the  row  algorithm  is  the  most  ef- 
ficient possible,  in  the  sense  that  there  exists  no  algorithm  to  construct 
S*  which  has  a smaller  order  of  computational  effort,  n(1og  n) , than  the 
row  algorithm.  For  example,  efficient  Implementations  of  the  first  and 
second  algorithms  In  [8]  result  respectively  In  computational  orders 
of  n^  and  n ; either  order  is  greater  than  that  of  the  row  al^otlthm.  For 
a discussion  of  a number  of  other  algorithms  of  order  n(log  n)  for  solving 
geometric  problems,  see  Shamos  [5  ],  [6  ]• 

Subsequently,  following  [8  ],  we  define  the  Line  Construction  Procedure: 

roughly  speaking  the  procedure  consists  of  plotting  the  points  P^ 

and  constructing  both  a vertical  and  horlsontal  line  through  each  point. 

If  the  algebraic  equivalent  of  the  line  construction  procedure  is  considered 
to  be  part  of  the  problem  formulation,  rather  than  part  of  the  row  algdrlthm, 
then  the  row  algorithm  is  of  order  n,  and  no  algorithm  to  construct  S*  can 
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be  of  an  order  smaller  than  n.  Even  when  the  line  construction  procedure  is 
considered  to  be  part  of  the  row  algorithm,  the  algorithm  performs  as  If  it 
is  of  order  n until  n becomes  "large." 

In  the  next  section,  after  presenting  and  Illustrating  a number  of 
definitions,  we  present  and  Illustrate  the  "arrow"  algorithm  which  constructs 
the  efficient  set.  We  then  give  a characterization  of  the  efficient  set. 

The  following  section  contains  the  row  algorithm.  The  last  section  of  the  paper 
consists  of  the  analysis  needed  to  justify  the  algorithms  and  the  characterization 


CHARACTERIZING  THE  EFFICIENT  SET 


Figure  1 Illustrates  a basic  notion,  a diamond  with  center 
and  radius  e^,  denoted  by  D(P^,  e^) . D(P^,  e^)  consists  of  all  points 

in  the  plane  whose  rectilinear  distance  from  P^^  is  no  greater  than  e^. 
The  boundary  of  D(P^,  e^^)  consists  of  all  points  in  the  plane  whose 
distance  is  equal  to  e^,  and  so  of  course  has  the  property  that  any  two 
points  on  the  boundary  are  the  same  rectilinear  distance  (e^)  from  P^. 

In  other  words,  the  boundary  is  a contour  line  of  the  rectilinear 
distance  from  P^  of  value  e^. 

We  call  line  segments  parallel  to  the  line  y ■ x SW  - NE  line 
segments , and  call  line  segments  parallel  to  the  line  y ■ -x  SE  - NW 
line  segments . Note  that  two  edges  of  any  D(P^ , e^)  are  SE  - NW  line 
segments,  while  the  other  two  edges  are  SW  - NE  line  segments. 

Line  Construction  Procedure.  Through  each  point  P^  construct  a 
horizontal  line  and  a vertical  line.  The  horizontal  (vertical), 
line  should  extend  at  least  as  far  right  and  as  far  left  (as  far  up 
and  as  far  down)  as  every  Pj^.  Subsequently  whenever  we  refer  to  a 
line  we  mean  a constructed  line  unless  we  specify  otherwise.  Figure  2a 
illustrates  the  construction  procedure. 

Noncollinearity  Assimption.  We  assume  that  not  all  the  P^  lie  on  a 
single  vertical  line,  or  on  a single  horizontal  line,  as  in  this  case 
S*  is  Just  the  line  segment  Joining  the  two  P^  which  are  farthest  apart, 
so  that  constructing  S*  is  a trivial  problem. 

Definitions.  Figure  3 illustrates  a number  of  the  definitions  to 


follow.  For  any  vertical  line  we  define  the  union  of  the  line  with 


the  set  of  points  t'  the  right  (left)  of  the  line  to  be  the  set  of  points 
which  are  east  (west)  of  the  line.  (Note  that  this  definition  permits 
a point  on  a vertical  line  to  be  both  east  and  west  of  the  line.)  Similarly 
we  define  the  set  of  points  north . and  the  set  of  points  south . of  each 
horizontal  line. 

Given  any  two  distinct  adjacent  horizontal  lines  H and  with  H 
north  of  H',  and  any  two  distinct  adjacent  vertical  lines  V and  V , with 
V east  of  V,  we  call  the  set  of  points  lying  west  of  V,  east  of  V,  south 
of  H,  and  north  of  H',  a box,  and  denote  the  box  by  B.  We  call  the  col- 
lection of  all  boxes  between  any  two  adjacent  vertical  (horizontal)  lines 
a column  (row) . Each  of  the  four  Intersections  of  the  box  with  a line 
we  call  an  edge  of  B.  We  say  two  boxes  are  adj acent  if  their  intersection 
is  an  edge  of  each  box.  The  collection  of  all  points  lying  south  of  H' 
and  east  of  V we  call  the  SE  direction  of  B (abbreviated  SE(B))  similarly 
we  define  and  directions  of  B,  and  use  the  abbreviations 

SW(B) , NW(B),  and  NE(B)  respectively.  We  call  the  abbreviations  SE,  SW, 

NW,  and  HE  the  box  direction  labels.  We  say  that  a direction  of  B Is 
unoccupied  (occupied)  If  there  is  no  (at  least  one)  In  the  direction. 

We  denote  the  union  of  all  the  boxes  by  B. 

Arrow  Drawing  Procedure.  For  each  box  B we  draw  an  arrow  pointing  from 
SE(B)  to  NW(B)  whenever  SE(B)  Is  unoccupied,  and  call  the  arrow  a SE 
arrow.  We  say  the  arrow  points  away  from  the  south  and  east  edges  of  B , 
and  points  towards  the  west  and  north  edges  of  B.  Likewise  we  construct 
and  define  SW,  NW,  and  NE  arrows  whenever  SW(B)  , NW(B),  and  NE(B) 
respectively  Is  unoccupied.  We  call  the  abbreviations  SE,  SW,  NW,  NE 
the  arrow  direction  labels.  Figure  2b  Illustrates  the  arrow  drawing 


procedure 


The  following  observation  provides  the  motivation  for  the  arrow 


drawing  procedure. 

Observation  1.  Suppose  we  are  given  any  box  B with  arrow  a In  B. 

For  any  point  X In  B,  move  the  point  X In  B In  the  direction  towards 
which  a points  along  a line  segment  parallel  to  a , until  the  point 
Intersects  a box  edge  at  a point  X’ . X and  X'  are  such  that 
r(X',  P^)  r(X,  ^ i ^ Further,  If  X'  X and  the  direction 

towards  which  a points  Is  occupied,  then  X'  dominates  X. 

To  establish  this  observation  (see  Figure  4 ),  suppoe?  without 

loss  of  generality,  that  a Is  a SE  arrow  and  X'  ^ X.  We  Viave 
r{X’,  P^)  - r(X,  P^)  for  P^  In  NE(B)  u SW(B),  while  r (X',  P^)  < r(X,  P^) 
If  P^  Is  In  NU(B).  Since  SE(B)  Is  unoccupied,  X'  thus  dominates  X 
whenever  NW(B)  Is  occupied. 

Observation  2.  Each  box  has  exactly  0,  1,  or  2 arrows.  Whenever  a 

box  has  2 arrows  the  arrows  are  parallel  and  point  In  opposite  directions. 

To  establish  this  observation  we  note  that  If  a box  has  perpendicu- 
lar arrows  then  there  are  no  P^  In  some  half-plane  defined  by  a line 
passing  through  some  edge  of  the  box,  which  Is  Impossible.  Hence  a 
box  has  at  most  2 arrows,  and  If  It  has  2 arrows  then  the  arrows  are 
parallel. 

Definitions.  We  call  a box  with  0,  1,  or  2 arrows  a null-box,  1-box, 
and  2-box  respectively.  For  any  1-box  B we  call  the  two  edges  towards 
which  the  arrow  In  B points  the  leading  edges  of  B (see  Figure  4). 

By  virtue  of  the  definitions  and  Observation  2,  we  have 
Observation  3.  Each  box  in  B Is  either  a null-box,  1-box,  or  2-box. 


A null-box  has  no  unoccupied  directions.  A 1-box  has  exactly  one 


unoccupied  direction,  which  has  the  same  direction  label  as  the  box 
arrow.  A.  2-box  has  exactly  two  unoccupied  directions,  which  differ 
by  180  degrees:  the  labels  of  the  unoccupied  directions  are  identical 
to  the  labels  of  the  two  arrows  in  the  box. 

Figure  5 illustrates  Observation  3. 

Observations  1 and  3 give 

Observation  4.  If  B is  a 1-box,  and  X is  any  point  in  B such  that  X 
is  not  on  a leading  edge  of  B,  then  there  is  a point  X*  on  a leading 
edge  of  B such  that  X'  dominates  X. 

It  is  also  convenient  to  state 

Observation  5.  Any  point  X which  is  not  in  B is  dominated  by  a point 
X’  which  lies  on  the  boundary  of  B,  and  is  the  closest  point  in  B to  X. 

We  now  state  the 

Arrow  Algorithm.  To  determine  the  set  S*  of  all  efficient  locations, 
carry  out  the  line  construction  and  arrow  drawing  procedures,  and  classify 
each  box  as  a null-box,  1-box,  or  2-box. 

If  there  are  no  1-boxes,  take  S*  “ B;  otherwise  choose  a 1-box  B 
not  yet  chosen  and  delete  from  B all  points  in  B except  those  on  the 
leading  edges  of  B;  repeat  this  deletion  procedure  for  every  1-box. 

Denote  by  B the  subset  of  B remaining  after  the  completion  of  the  de- 
letion procedure.  Take  S*  ■ 6. 

We  remark  that  if  E is  a common  edge  of  two  1-boxes , B and  B ' , 
if  E is  a leading  edge  of  B,  and  not  a leading  edge  of  B',  then  E (except 
for  one  endpoint)  will  be  deleted  from  B once  B'  is  chosen.  Figures 

2a,  2b,  and  2c  Illustrates  the  algorithm:  a null-box  is  identified  by  a 
dot  in  the  box. 
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Due  to  Observations  4 and  5,  we  have 


Observation  6.  Any  point  not  In  B la  dominated  by  a point  In  B. 

With  two  additional  definitions  we  can  characterize  S*. 

Definition.  An  edge  E of  any  box  B Is  called  a connecting  edge  If 

(a)  every  arrow  In  B points  towards  E and  E Is  contained  In  the 
boundary  of  B,  or 

(b)  there  Is  also  a box,  say  B',  such  that  E Is  the  common 
edge  of  B and  B',  and  every  arrow  In  each  box  points  towards 
E. 

I 

Connecting  edges  are  Illustrated  In  Figure  6.  | 

Definition.  Denote  by  6*  the  union  of  the  following; 

(a)  all  null-boxes 

(b)  all  2-boxes 

(c)  all  connecting  edges. 

Subsequently  we  establish  that  S*,  the  set  of  efficient  points, 

B,  the  set  of  points  left  by  the  algorithm,  and  B*  are  all  Identical. 

Thus  the  arrow  algorithm  deletes  all  points  In  B which  are  not  points 

In  B*.  Figure  2c  Illustrates  B*.  ^ 

As  a final  comment,  we  note  that  the  arrow  algorithm  as  stated,  | 

and  each  of  the  algorithms  in  [8  ].  is  "memoryless"  to  some  extent,  | 

In  the  sense  that  each  algorithm  can  make  more  use  than  It  does  of 
Information  obtained  during  the  process  of  determining  the  efficient 
set.  In  the  following  section  we  present  the  row  algorithm,  which  ex- 
ploits such  Information  efficiently. 


CONSTRUCTING  THE  EFFICIENT  SET: 

THE  RON  ALGORITHM 

Some  notation  Is  convenient.  Denote  the  horizontal  lines  by  H^,  H^. 

H from  north  to  south,  and  the  rows  by  R, , ...»  R from  north  to  south, 
p+l  1 p 

9 

For  1 1 ^p+1,  denote  by  x coordinate  of  the  westmost  (eastmost) 

existing  facility  location  on  H^.  For  R^,  1 ^ 1 ^P»  define  NW^(SW^)  to  be 
the  X coordinate  of  the  westmost  existing  facility  location  which  Is  north 
(south)  of  R^.  Define  NE^(SE^)  to  be  the  x coordinate  of  the  eastmost 
existing  facility  location  which  is  north  (south)  of  R^.  (See  Figure  11.) 

As  an  Immediate  consequence  of  the  definitions  we  have 
Observation  7 Let  B be  a box  in  R^^,  1 1 ^ p* 

(a)  NW(b)  [SN(b)]  Is  unoccupied  If  and  only  If  B Is  west  of  the  vertical 
line  X = NW^(x  ■ SW^^) . 

(b)  NE(B)  [SE(B)]  Is  unoccupied  If  and  only  If  B Is  east  of  the  vertical 
line  X - NE^(x  - SE^) . 

By  virtue  of  the  above  observation  we  can  readily  classify  the  boxes 
in  each  row  R^.  The  following  observation  facilitates  the  computations  of 
NWi,  SW^,  NE^,  and  SE^^. 

Observation  8 The  following  recursive  relationships  are  true: 


NWj^ 

- w 

1 

NEj^ 

-=1 

NW^ 

- mln(W^,  NW^_^) 

NEi 

- max(E^,  t 2 ^ 1 ^ p 

SWi 

- mln(W^j.  SW^^j^) 

SE^ 

- im:i(E^^j^,  SE^_j^j^)  , 1 <_  1 <_  p-1 

sw 

■ M 

SE 

■ E , 

p 

p+l 

P 

p+l. 

9 


The  above  recursions  are  easily  established.  For  example,  certainly 
NWj^  “ Also,  the  only  existing  facility  locations  north  of  lie  on 

^1’  ^2’  ^i’ 

■ mln(Mj^,  W^,  ...,  W^) . 

Likewise 

NW^_j  • mln(Wj^,  '^i-l^ 

and  so  certainly 

NW^  - mln(W^,  NW^_j^) . 

An  equivalent  geometric  means  of  computing  the  recursions,  which 

both  provides  insight  and  is  easy  to  carry  out  manually,  may  be 

described  as  follows. 

\ 

4-Color  Procedure.  Associate  the  colors  Blue,  Green,  Red,  and  Yellow 
with  the  directions  NW,  NE,  SW,  and  SE  respectively. 

Repeat  the  follwolng  North  to  South  Step  for  lines  ...,  H , consecutively. 

p+1 

Beginning  at  the  west  (east)  boundary  of  S,  draw  a blue  (green)  line  over 

from  west  to  east  (east  to  west)  terminating  the  line  at  x ■ (x  ■ NE^) , which 

is  the  point  at  which  one  of  the  following  two  events  first  occurs: 

(a)  The  line  initially  Intersects  an  existing  facility  location  on  H^, 

(b)  The  line  attains  the  same  length  as  the  blue  (green)  line  on 

Repeat  the  following  South  to  North  Step  for  lines  ...,  Hi  consecutively. 

Beginning  at  the  west  (east)  boundary  of  B,  draw  a red  (yellow)  line  over 
from  west  to  east  (east  to  west),  terminating  the  line  at  x ■ SW^  (x  ■ SE^^)  , which 
is  the  point  at  which  one  of  the  following  two  events  first  occurs : 

(a)  The  line  initially  Intersects  an  existing  facility  location  on  H^, 

(b)  The  line  attains  the  same  length  as  the  red  (yellow)  line  on 

Note,  in  the  north  to  south  step  that  when  an  existing  facility  lying  on 
the  west  (east)  boundary  of  0 and  is  first  encountered,  every  blue  (green) 
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line  drawn  subsequently  is  a degenerate  line  of  zero  length.  In  the 
south  to  north  step  when  an  existing  facility  lying  on  the  west  (east) 
boundary  of  B and  Is  first  encountered,  every  subsequent  red  (yellow) 
line  dram  is  a degenerate  line  of  zero  length.  Figure  11  illustrates  the 
use  of  the  4-color  procedure. 

Classifying  the  Boxes. 

Once  the  4-color  procedure  has  been  completed,  if  we  suppose  the  initial 
horizontal  lines  to  be  uncolored,  it  is  easily  verified  that  the  north  edge 
of  each  box  colored  in  the  north  to  south  step  is  exactly  one  of  the  colors 
blue  or  green,  while  the  south  edge  of  each  box  colored  in  the  south  to  north 
step  is  exactly  one  of  the  colors  red  or  yellow.  A box  is  a NW  1-box  (NE  1-box) 
if  and  only  if  its  north  edge  is  blue  (green)  and  its  south  edge  is  either  un- 
colored or  blue  (uncolored  or  green) . A box  is  a SW  1-box  (SE  1-box)  if  and 
only  if  its  south  edge  is  red  (yellow)  and  its  north  edge  is  either  uncolored 
or  red  (uncolored  or  yellow) . A box  is  a NW-SE  (NE-SW)  2-box  if  and  only  if 
its  north  and  south  edges  are  blue  and  yellow  (green  and  red)  respectively. 

A box  is  a null-box  if  and  only  if  both  its  north  and  south  edges  are  un- 
colored. The  color  combinations  listed  in  this  paragraph  for  the  three  types 
of  boxes  are  the  only  ones  possible. 

Definitions. 

(1)  We  observe  that  NW^(SW^)  is  the  projection  onto  the  x axis  of  the  east 
tip  of  the  blue  (red)  line  on  ^E^  (SE^)  is  the  projection 

onto  the  x axis  of  the  west  tip  of  the  green  (yellow)  line  on 
refer  to  NW^,  SW^,  NE^,  and  SE^  as  the  blue,  red,  green,  and  yellow  projections 
for  row  i.  For  each  row  i,  we  define 
W*  - max  (NW^,  SW^) 

E*  - min  (NE^,  gg^) 

and  call  H*  and  E*  the  west  and  east  projections  for  row  1 . See  Figure  11  for 
an  example. 


(ii)  Given  a box  B and  distinct  vertical  lines  x * V'  and  x ~ V with  V'  west 
of  V,  we  write  V < B < V when  we  mean  V'  is  west  of  B and  B is  west  of  V. 


Note  tliat  for  a given  row,  there  is  at  least  one  box  B in  the  row  such  that 
V'  < B < V if  and  only  if  V'  < V.  Further,  for  a given  box  B and  vertical  line 
V",  either  B < V"  or  V"  < B. 

* it 

(Hi)  Comparisons  of  and  for  are  crucial  to  the  row  algorithm  we 
shall  develop  for  finding  S*.  Define  to  be  in  Condition  0 (C-0),  Con- 
dition 1 (C-1),  or  Condition  2 (C-2)  when  W*  < E*,  W*  - E*,  and  W*  > E* 
respectively.  Note  the  three  conditions  are  mutually  exclusive  and  ex- 
haustive. 

(iv)  For  each  line  H^,  1 £ i £ p + 1,  define 

u^  « max(NW^,  SW^_^)  v^  - min(NE^,  SE^_^) 

where,  by  convention,  SW  ■ - » SE  ■ “>. 

0 0 

Interpreting  the  Conditions. 

(C-0)  Since  a box  in  R^  is  a null-box  if  and  only  if  its  north  and  south 
edges  are  uncolored,  the  following  conditions  are  all  equivalent  for  at 

least  one  box  B in  R^  to  be  a null-box:  the  projection  of  B lies  between 

* * * * 

the  west  and  east  projections  for  R^j^;  < B < E^;  < E^^;  R^^  is  in  C-0. 

* * a * 

When  < E^^  the  null-boxes  in  R^^  are  those  boxes  B for  which  < B < E^^. 

(C-2) . Since  a box  in  R^  is  a 2-box  if  and  only  if  Its  north  and  south  edges  are 
either  blue  and  yellow  or  green  and  red  respectively,  there  is  a least  one 

2-box  in  R^  if  and  only  if  the  west  and  east  projections  for  R^  overlap,  that 

* * 

is,  > E^,  that  is,  R^  is  in  C-2.  Subsequently  we  establish  there  is  ex- 
actly one  2-box  B in  R^  if  and  only  if  E*  < B < U*. 

* * 

(C-1) . Denote  by  V the  vertical  line  identical  to  the  line  x ■ > E^.  Denote  by 

E that  part  of  V lying  south  of  and  north  of  Since  W*  ■ E*  the  west 

and  east  projections  meet,  so  each  of  the  boxes  of  which  E Is  a vertical  edge 
has  at  least  one  colored  edge.  If  E is  an  oast  vertical  edge  of  B,  at  least 
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Row  Algorithm. 

(1)  Rank  the  existing  facilities  by  their  y coordinates  to  determine  the 

lines  H, , ....  H . , . 

1 p+l 

(2)  Compute  and  l^l^p  + l. 

(3)  Compute  NW^  and  NE^  using  the  recursions  of  Observation  8,  1 i £ P* 

(4)  Compute  SW^  and  SE^  using  the  recursions  of  Observation  8,  1 f.  1 P* 

* * 

(5)  Compute  E^^  and  for  1 1 P* 

. * * 

(6)  (a)  < E^  If  and  only  If  Is  In  C-0.  The  collection  of  null-boxes 

In  R^  Is  the  collection  of  boxes  B In  R^^  for  which  W*  < B < E*. 

Place  all  null-boxes  In  B*. 

* * 

(b)  ■ E^  if  and  only  If  Is  In  C-1.  The  unique  vertical  connecting 

* * 


edge  E In  R^  Is  the  vertical  edge  E contained  In  the  line  x ■ > E^. 

Place  E in  3*. 

* * 


(c)  Wj  > Ej^  If  and  only  if  R^  Is  in  C-2.  The  unique  2-box  in  R^^  Is  the 
box  B such  that  E*  < B < W*;  place  this  box  in  3*.  All  null-boxes,  2 
boxes,  and  vertical  connecting  edges  are  now  In 

(7)  For  the  horizontal  line  H^,  1 f.  1 £ p + 1; 

(a)  “ Vj  then  contains  no  horizontal  connecting  edge; 

(b)  If  ^ ® horizontal  edge  E contained  In  Is  a horizontal 


connecting  edge  if  and  only  if  it  lies  in  between  the  vertical  lines 


X ■ u^  and  x ■ v^,  and  Is  not  an  edge  of  any  null-box.  Place  each  such  E 

(8)  S*  - 3*.  Stop. 

Figure  11  illustrates  the  use  of  steps  2 through  6 of  the  algorithm, 
together  with  the  result  of  step  7. 

CosB>utatlonal  Efficiency  Questions. 
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* 
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the  base  2,  and  n is  the  ntinber  of  existing  facilities.  We  shall  see  that 
the  order  of  computational  effort  of  the  row  algorithm  is  n(log  n) . Hence 
the  row  algorithm  is  "optimal",  in  the  sense  that  its  order  of  computational 
effort  is  as  small  as  that  of  any  algorithm  for  constructing  S*. 

Given  n distinct  nundiers,  a^^,  ...»  a^,  it  is  known  that  (see  pps.  159-170 
of  [3],  or  pps.  65-67  of  [1]),  the  minimum  number  of  comparisons  needed  to  rank 
the  numbers  is  of  order  n(log  n) . Now  define  existing  facility  locations 
Pj  “ (Sj,  Sj),  1 f.  j £ n,  and  consider  the  example  problem  of  specifying  S* 
for  the  existing  facility  locations.  In  order  to  state  S<*,  let  [1],  ...» 

[n]  be  a permutation  of  1,  2,  ...,  n such  that 
*(J1  ^ J define 

S*-((x.  y),  ajj,  < X.  y < liJ  <n-  1. 

* 

Each  Sj  is  a 2-box  with  SW  and  NE  vertices  of  and  respectively. 

It  is  easy  to  use  the  row  algorithm  to  establish  Chat 

* * * 

S*  • S,  u S-  u . . . u S , . 

1 L n— i 

Any  other  correct  algorithm  to  determine  S*  would  also  need  to  specify 

* * 

each  Sj , and  each  Sj  ia  completely  determined  by  the  nundiers  a j j j and 

*[j+ll’  ®Ij]  ^ *Ij+l]*  ^ ^ **ence  in  order  to  specify  S*  any 

correct  algorithm  must  be  able  to  sort  the  numbers  a^ a^,  and 

so  its  order  of  computational  effort  is  at  least  n(log  n) . (Mote  this 

example  problem  also  establishes  that  any  algorithm  to  construct  S*  is 

at  least  of  order  n when  the  Line  Construction  Procedure  is  considered  not 

to  be  part  of  the  algorithm, as  S*  is  the  union  of  n - 1 boxes.) 

Now  consider  Che  row  algorithm.  Step  (1)  requires  the  ranking  of  n 

numbers,  and  so  is  of  order  n(log  n) . If  n^  is  the  number  of  existing 

facilities  on  H^,  finding  W^  and  is  equivalent  to  the  problem  of  find- 
ing the  maximim  and  miniisim  numbers  in  a collection  of  n^  nimd>ers,  and  it 
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Is  known  (see  Pohl  [4])  that  this  problem  can  be  solved  in  <(1.5)nj^  - 2> 
comparisons  at  best  (here  <y>  denotes  the  smallest  Integer  no  leas  than 
y) . Since  <(1.5)nj^  ” 2>  (1.5)  (n^^  - 1),  an  upper  bound  on  the  number  of 

comparisons  made  In  Step  (2)  Is  given  by 

[ {(1.5)(n^  - 1):  1 < 1 < p + 1}  - (1.5)n  - (1.5)(p  + 1). 

Steps  (3)  through  (6)  each  require  2p  comparisons,  for  a total  of  8 p com- 
parisons. In  the  analysis  (Observation  20)  we  establish  that  Step  (7)  can 
be  accomplished  in  no  more  than  8(p  + 1)  comparisons.  Thus  Steps  (2) 
through  (7)  require  at  most  (1.5)  n - (1.5) (p  + 1)  + 8p  + 8(p  + 1)  - 
(1.5)  n -f  (6.5)  (p  -f  1)  + 8p  £ 16n  comparisons.  The  minimum  number  of  com- 
parisons needed  for  Step  (1)  is  n <log  n>  - 1,  which  Is  no  greater 

than  n <log  n>.  Thus  an  upper  bound  on  the  number  of  comparisons  made  In 

Steps  (1)  through  (7)  Is  given  by  n <log  n>  + 16n.  Therefore,  roughly 
speaking.  Step  (1)  requires  more  effort  than  Steps  (2)  through  (7)  when 
n <log  n>  £ 16n,  which  Is  true  when  n £ 2^^  "65 ,536.  Hence  the  nonlinear 
effort  of  the  first  step  does  not  predominate  until  n Is  "large".  In  which 
case  the  row  algorithm  is  of  order  n(log  n) . In  effect,  setting  up  the  pro- 
blem (Step  (1))  requires  more  effort  than  solving  it  when  n Is  large! 

Some  Idea  of  how  the  computational  effort  of  the  row  algorithm  compares 

2 3 

with  algorithms  of  order  n or  n can  be  obtained  by  examining  Table  1.  The 

2 3 

comparisons  In  the  table  are  conservative,  in  the  sense  that  order  n and  n 
algorithm  typically  have  a multiplicative  constant  greater  than  one. 
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ANALYSIS 


Diamond  Intersection  Lemma.  Given  diamonds  D(Pj^,  i ® I ^ ^ 

with 

D=  n {D(P^,  e^)  : 1 e 1}  <|> 

If  X Is  a point  In  D such  that  no  point  in  D dominates  X,  than  X is  in  S*. 
Proof.  Let  Y be  any  point  in  the  plane.  If  Y / D there  is  an  index  q e I 
such  that  Y / D(P„,  e^),  and  thus  r(P„,  X)  < e < r(P  , Y),  so  Y does  not 

q q q — q q 

dominate  X.  By  hypothesis  Y does  not  dominate  X if  Y e D,  so  X e S*.  (We 
shall  refer  to  this  lemma  as  DIL.) 

Property  1.  Every  null-box  is  efficient. 

Proof.  Let  B be  any  null-box  and  let  X e B.  Since  B Is  a null-box  we 

may  choose  P^^,  P2.  P^  and  P^  in  NE(B) , NW(B),  SW(B) , and  SE(B)  respectively 

as  Illustrated  In  Figure  7 . Define  e^  = r(X,  P^) , and  construct 

Di  E D(Pj^,  e^),  1 1 1 £ 4.  With  Lj^^  E D^^  n D^,  L^^  = D2  n D^,  we  note 

L, . is  a SE-NW  line  segment,  L Is  a SW-NE  line  segment,  and 
IJ  24 

D - L^3  n L2^  - {X}. 

Thus  no  point  In  D dominates  X.  DIL  thus  implies  X e S*.  Since 
X Is  an  arbitrary  point  In  B,  B c s*. 

Observation  9.  Only  a 1-box  can  be  adjacent  to  a 2-box.  Whenever  a 
1-box  and  2-box  are  adjacent,  the  arrows  In  the  boxes  are  parallel, 
and  the  arrow  In  the  1-box  points  towards  the  edge  comnon  to  the  1-box 
and  2-box. 

Proof.  Let  B be  a 2-box,  and  let  B'  be  a box  adjacent  to  B.  It  la 
direct  to  verify  that  If  B'  is  a null-box,  an  occupied  direction  of 
B'  is  contained  In  an  unoccupied  direction  of  B,  which  la  Impossible. 
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If  B*  Is  a 2-box,  either  an  occupied  direction  of  one  box  Is  contained  In 


an  unoccupied  direction  of  the  other,  or  else  there  Is  no  on  the  line 
separating  the  boxes;  both  Impossible  situations.  If  B*  is  a 1-box,  an 
occupied  direction  of  B'  is  contained  in  an  unoccupied  direction  of  B ex- 
cept when  the  arrows  In  the  two  boxes  are  parallel,  which  completes  the 
proof. 

Observation  10.  Let  B be  a 2-box,  and  let  B'  be  a 1-box  adjacent  to  B,  and 
in  the  same  row  (column)  of  boxes  as  B.  If  B"  is  any  other  box  in  the 
same  row  (column)  as  B and  B*  such  that  B'  lies  between  B and  B",  then  B" 
is  also  a 1-box,  and  the  arrows  in  B*  and  B"  have  the  same  label. 

Proof.  Let  B be  a 2-box.  By  Observation  9,  we  know  that  any  box  B'  ad- 
jacent to  B is  a 1-box.  We  may  assume  that  B,  B',  and  B"  are  all  in  the 

same  row.  Also,  due  to  symmetry,  we  may  assume  that  the  arrow  labels  of  B 
are  NE  and  SW,  and  that  B'  is  east  of  B.  Thus  B"  will  be  east  of  B’,  as  Illustra- 
ted in  Figure  8.  Now  SW(B')  c SW(B")  and  SW(B')  occupied  Implies  SW(B")  Is  oc- 

cupied. NW{B)  c NW(B")  and  NW(B)  occupied  implies  NW(B")  is  occupied. 

NE(B")  c NE(B)  and  NE(B)  unoccupied  Implies  NE(B")  is  tmoccupied.  If  SE(B'') 
is  unoccupied,  no  P^  would  lie  on  the  vertical  line  passing  through  the 
east  edge  of  B",  which  is  Impossible.  Thus  SE(B")  is  occupied,  and  B"  is 
a 1-box  having  an  arrow  with  a NE  label,  the  same  label  as  the  arrow  In  B* . 
Property  2.  Every  2-box  is  efficient. 

Proof.  Let  B be  any  2-box.  Without  loss  of  generality  we  may  assume  the 
arrows  in  B are  NE  and  SW  arrows,  that  B is  the  intersection  of  half-planes 
defined  by  lines  as  illustrated  in  Figure  9.  If  there  is  a 1-box  B*  which 
is  adjacent  and  north  of  B,  by  Observations  9 and  10,  the  arrow  in  B' 
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has  a NE  label.  If  there  Is  an  adjacent  box  B"  to  B such  that  B"  is  south 
of  B,  then  the  arrow  in  B"  has  a SW  label.  Since  SW(B)  is  unoccupied,  there 
roust  be  some  e V n NW(B) , or  else  there  would  be  no  e V.  Likewise, 
there  must  be  some  P^  e V n SE(B) . Choose  any  X e B,  and  define 

e^  - r(X,P^),  ej  - r(X,Pj),  D - D(Pj , e^)  n D(Pj , e^).  Since  X e D, 

by  OIL  X e S*  provided  no  point  in  D dominates  X.  We  make  the 
following  observation:  D is  a SW-NE  line  segment,  entirely  contained 
in  the  column  of  boxes,  say  C,  which  contains  B and  B'.  Now  if  Y e D n B, 
r(X,  Pj^)  ■ r(Y,  Pj^)  for  all  i,  so  Y does  not  dominate  X,  thus  completing 
the  proof  for  the  case  when  B'  * ^ ■ B".  If  Y e D,  Y / B,  Y is  in  a box 
in  C which  is  a 1-box,  and,  hence,  by  Observations  9 and  10,  there  is  a 
point  Y'  in  B n D which  dominates  Y.  Since  Y'  dominates  Y,  and  Y'  does 

not  dominate  X,  Y does  not  dominate  X,  so  X e S*.  Thus  B c S*. 

Property  3 Every  connecting  edge  is  efficient. 

Proof . Let  E be  any  connecting  edge.  Without  loss  of  generality,  we  may 
assume  E is  an  east  edge  of  a 1-box  B,  and  that  B has  a SW  arrow. 

Consider  first  the  case  where  E is  contained  in  the  boundary  of  6. 

Then,  as  Figure  6a  illustrates,  since  E lies  on  the  east  boundary  of  B, 

SW(B)  is  unoccupied,  and  there  is  a point  e H n V.  Also,  since  E lies 

on  the  east  boundary  of  B and  NE(B)  is  occupied,  there  is  a point  e 
V n NE(B).  If  X E E,  e^  = r(X,  P^^),  e^  = r(X,  Pj^),  and  D = D(P^,  e^)  n 
D(Pk>  then  D ■ {X},  and  so  by  DIL,  X e S*  and  thus  E c S*. 

Now  suppose  there  is  a 1-box  B',  adjacent  to,  and  east  of,  B.  The 
arrow  in  B'  is  either  a SE  arrow,  as  illustrated  in  Figure  6b,  or  a NE 
arrow,  as  Illustrated  in  Figure  6c,  in  the  torm»r  case  we  conclude  there 
is  some  P^  E SE(B)  n SW(B*),  some  P^  e V n NW(B),  end  some  Pj^  e V”  n NE(B'), 
as  illustrsted  by  Figure  6b.  Hence  if  X e E,  defining  e^  - r(X,  P^), 
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X must  lie  In  one  such  box.  If  the  box  is  a null-box  or  2-box,  certainly 
X ( 3*,  so  It  remains  to  consider  the  case  where  X is  In  a 1-box,  but  X is 
in  no  null-box  or  2-box. 

Since  X is  in  a 1-box,  say  B^,  and  in  B,  X must  lie  on  some  leading 
edge  of  B^,  say  E,  of  positive  length.  If  E is  contained  in  the  boundary 
of  3,  since  the  arrow  in  points  towards  E,  E is  a connecting  edge.  Thus 
it  remains  to  consider  the  case  where  E is  not  contained  in  the  boundary  of 
3.  In  this  case  there  is  some  box,  say  b^,  such  that  E - Bj^  n B2.  Thus 
X is  in  B2,  and  so  B2  must  be  a 1-box.  If  X is  not  an  endpoint  of  E,  since 
X e 8 the  arrows  in  B^  and  B2  must  point  towards  E,  and  so  E is  a connecting 
edge.  Thus  suppose  X is  at  an  endpoint  of  E.  Since  X e B^  n B2  and  X is 
not  on  the  boundary  of  3,  there  exist  1-boxes  B^  and  B^  such  that  {X}  ■ 

B^  n B2  n B^  n B^  (as  illustrated  in  Figure  10  ).  Thus,  by  Observation  11, 

X lies  on  a connecting  edge. 

Theorem  1.  3 . 3^,  and  S*  are  identical  and  nonempty. 

Proof . By  Observation  6 , every  point  not  in  8 is  dominated,  and  thus 
S*  c 3.  By  Properties  1,  2,  and  3,  8*  c S*,  and  so  8*  c S*  c 8 . Since 
Property  A gives  8 =8**  we  conclude  8*  ■ S*  ■ 6.  As  every  P^^  is  in  S*, 

S*  ^ A,  so  8*  - 8 - S*  A. 
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We  now  give  the  analysis  needed  to  Justify  the  row  algorithm. 

* * 

Observation  12.  A box  B in  is  a null-box  if  and  only  if  < B < 

Proof.  By  Observation  7«  the  four  directions  of  a null-box  B are  oc- 
cupied if  and  only  if  < B,  < B,  B < NE^^,  B < SE^,  which  is  equlv- 

* * 

alent  to  < B < E^. 

Observation  13.  A box  B in  is  a 2-box  if  and  only  if  either  NEj^  < B < SE^ 
or  NWj  < B < 

Proof.  By  Observation  3,  a box  B in  R^  is  a 2-box  if  and  only  if  B is 

either  a NW-SE  or  a NE-SW  2-box.  B in  is  a NE-SW  2-box  if  and  only  if 

NE(B)  and  SW(B)  are  unoccupied  and  NW(B)  and  SE(B)  are  occupied.  Thus 
Observation  7 implies  B is  a NE-SW  2-box  if  and  only  if  NW^  < B,  NE^^  < B, 

B < SW^,  B < SE^.  By  definition,  NW^^  £ NE^  and  SW^^  £ SE^^,  so  B in  R^  is 
a NE-SW  2-box  if  and  only  if  NE^^  < B < SW^^.  Similarly,  B in  Rj^  is  a NW-SE 
2-box  if  and  only  if  SE^^  < B < NW^^. 

Observation  14.  Given  a box  B in  R^, 


NEj  < B<  SWj^ 

(1) 

or 

SE^  < B < NW^ 

(11) 

If  and 

only  If 

* * 

(111) 

Ef  < B < Wf  . 

Proof. 

By  definition,  E*  £ SE^,  E*  £ NE^,  SW^  £ W*,  and  NW^^  £ W*, 

so  If 

(i)  or  (li)  is  true  then  (Hi)  is  true. 

Suppose  (ill)  is  true.  Either  E^  ••  SE^  or  ■ NE^.  When  E^  ■ SE^, 
(ill)  and  SW^  < SE^  give 

SWi  < SEi  - eJ  < B < max(SW^.  NW^) . (Iv) 

Since  (iv)  implies  B < NW^,  we  conclude  (111)  Implies  (11).  When  E^  ■ NE^, 
(111)  and  NW^  < NE^  give 


* * 

the  fact  that  E is  a vertical  connecting  edge.  Thus  ■ E . 

h h a * 

It  reaains  to  show  - V - E^.  Let  x - V’  be  a line  such  that  - V - E^ 

and  suppose  V ^ V.  Since  W*  ■ V , every  box  west  of  V'  will  have  either  a 
NW  or  SW  label,  while  since  V'  ••  E^  every  box  east  of  V'  will  have  either  a 
NE  or  SE  label.  Denote  by  E'  the  vertical  edge  in  contained  in  V . From  the 
first  part  of  the  proof  we  know  E’  is  a vertical  connecting  edge.  Without  loss 
of  generality  we  may  assume  V'  Is  west  of  V.  Thus  there  exists  a 1-box  B'  of 
such  that  E'  Is  a west  edge  of  B' , and  there  exists  a 1-box  B of  R^  such  that 
E is  an  east  edge  of  B.  Since  E and  E'  are  vertical  connecting  edges,  there  Is 
a NW  or  SW  arrow  In  B pointing  towards  E,  and  a NE  or  SE  arrow  In  B'  pointing 
towards  E'.  Further,  NE(B)  and  SE(B)  are  occupied  while  either  SE(B')  or  NE(B')  Is 
unoccupied.  But  NE(B)  and  SE(B)  are  contained  in  SE(B')  and  NE(B')  respectively. 
Thus  an  occupied  direction  of  B is  contained  in  an  unoccupied  direction  of 
B'  if  V f*  V’ , and  we  have  a contradiction.  Thus  W*  ■ V ■ E*. 

Given  any  horizontal  edge  E,  whenever  E Is  east  of  the  vertical  line 
x - V and  west  of  the  vertical  line  x - V'  we  write  V < E < V . Given  any 
edge  E contained  In  It  Is  easy  to  verify  that  E Is  uncolored  If  and  only  If 
u^  < E < v^.  Hence  the  following  result  can  be  established: 

Observation  18.  Given  an  edge  E contained  In  a horizontal  line  H^,  E Is  a 
horizontal  connecting  edge  If  and  only  If  u^^  < E < v^^  and  E la  not  an  edge  of 
any  null-box. 

We  now  state 

Theorem  2.  The  set  constructed  by  the  row  algorithm  la  the  efficient  set. 

Proof.  A direct  consequence  of  Observations  12,  IS,  17,  18,  and  Theorem  1. 

It  now  only  remains  to  consider  the  computational  effort  involved  in  Step  (7) 
of  the  row  algorithm.  Due  to  Observation  18,  if  u^  ~ v^  there  is  no  connecting 
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edge  in  H . If  u.  < v.  and  we  remove  from  H.  the  interior  points  of  all 


edges  which  are  null-box  edges,  then  the  remaining  edges  (if  any)  are  the 


horizontal  connecting  edges  in  H. . We  consider  the  case  where  there  are 


null-boxes  in  both  ^ and 


as  this  case  requires  the  most  effort.  Let 


E. ].  I is  the  X projection  of  the 


uncolored  portion  of  H . , while,  by  Observation  12,  J'  and  J"  are  the  x-projections 


Because  the  null-boxes  whose 


horizontal  edges  are  in  H.  have  these  edges  uncolored,  we  have  J'  c I,  J"  c I 


The  most  direct  way  to  determine  the  computational  effort  for  this  case  is 


to  state  a simple  algorithm  which  identifies  all  the  horizontal  connecting  edges 


in  H . We  consider  the  algorithm  self-evident,  as  it  simply  determines  all  the 


points  in  I which  are  not  interior  to  J*  u J' 


Observation  19.  When  R 


and  R.  each  contain  null-boxes,  the  horizontal  con- 


necting edges  in  H.  may  be  determined  as  follows 


(a)  check  to  see  tdiich  of  the  intervals  J’  and  J"  has  a leftmost  endpoint; 


denote  this  interval  by  I'  and  denote  the  other  interval  by  I".  Let  1 


(b)  Check  to  see  if  I'  abuts  the  left  endpoint  of  I;  if  it  does  not,  the  edges 


in  H whose  x projections  11«  in  [u 


a' ] are  horizontal  connecting  edges 


(c)  Check  to  see  if  I'  and  I"  intersect 


X projections  lie  in  [b 


a"]  are  horizontal  connecting  edges 


(d)  Check  to  see  which  of  the  two  intervals  I'  and  I"  has  a rightmost  right 


endpoint,  and  denote  this  interval  by  1 


(e)  Check  to  see  if  I*"  abuts  the  right  endpoint  of  1;  if  it  does  not,  the  edges 
in  whose  x projections  lie  in  [b**,  v^]  are  horizontal  connecting  edges. 

In  steps  (a)  through  (e)  respectively  we  note  that  the  following  teme  are 


I 


compared:  W and  W u.  and  a';  b'  and  a”;  b'  and  b";  b"'  and  v..  Hence, 

1“  A 1 ^ ^ ^ 

given  all  the  necessary  data,  we  can  find  all  the  horizontal  connecting  edges 

by  making  at  most  5(p  + 1)  comparisons.  To  obtain  the  data,  with  reference 

to  Step  (7),  we  must  compute  the  u^  and  (2  p comparisons)  and  then  compare 

them  (p  + 1 comparisons).  Thus  Step  (7)  requires  at  most 

5(p  + 1)  + 2 p + (p  + 1)  ^ 8(p  + 1)  comparisons.  As  p + 1 £ n,  we  have 

Observation  20.  Step  (7)  of  the  row  algorithm  requires  at  most  8(p  4-  1)  j<  8 n 

comparisons. 


n 

2 

n 

3 

n 

r(n) 

2 

A 

8 

33 

A 

16 

6A 

69 

8 

6A 

512 

1A5 

16 

256 

A,  096 

305 

32 

1,02A 

32,768 

6A1 

6A 

A,  096 

262, lAA 

1,3A5 

128 

16,38A 

2,097,152 

2,817 

256 

65,536 

16,777,216 

5,889 

Table  1:  Computational  Effort  Comparisons  for  n , n , and  Row  Algorithm, 


r(n)  = n < log  n > - 


2<  log  n > 


+ 1 + 16  n 
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